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1.0  INTRODUCTION 


Elemental  boron  has  long  been  ol  interest  as  an  advanced  propellant  because  ot  its  high 
energy  density.  Recently,  there  has  been  an  increased  interest  in  paniculate  boron  oxidation  in 
tluonne  eniiched  environments.  This  is  due.  in  part,  to  applications  ot  dinuoroamino/niiroamino 
ba^ed  oxidizers  and  new  difluoroamino/azido  oxeiane  binders  and  solid  boron  in  advanced 
underwater  explosives,  new  propellant  tormulations  and  solid  propellant  igniuon  systems.  The 
potential  advantages  ot  fluorine  oxidation  include  elimination  ot  condensed  phase  boron  oxides  and 
oxyhydrides.  as  well  as  taster  bum  times.  Additionally,  fluorine  has  been  tound  to  enhance 
boron's  ignition  characteristics  in  oxygen  environments  by  promoting  gasification  of  the 
ubiquitous  oxide  layer. 

Considering  molecular  oxygen  and  molecular  fluorine  as  oxidizers,  it  is  also  interesting  to 
note  (see  Table  1)  that  the  heat  of  reaction  in  forming  gaseous  BF3  is  higher  or  comparable  to  the 
energy  released  in  oxidation  to  liquid  B2O3.  The  latter  is  difficult  to  achieve  within  typical 
residence  times  for  practical  devices  because  of  a  kinetic  bottleneck  in  the  condensation  process  of 
B203(g)  to  B203(1)  through  formation  of  the  intermediate  HB02(g).'  “  In  mixtures  con.sisting  of 
both  fluorine  and  oxygen,  the  equilibrium  boron  product  distribution  shifts  toward  OBF. 
Energetically,  OBF  formation  is  comparable  to  B203(l)  based  on  a  mass  weighting  of  the  fuel 
only;  however,  its  formation  is  less  desirable  when  based  on  the  total  weights  of  fuel  and  oxidizer. 
Heats  of  reaction  will  vary  with  different  fluorinated/oxygenated  oxidizers.  Because  the  heat  of 
formation  for  HF  is  more  exothermic  than  that  for  H2O,  the  overall  heal  released  from  combined 
B/O/H/C/F  based  systems  will  generally  be  competitive  with  B/O/H/C  based  systems. 


Table  1 .  Heats  of  Oxidation  and  Fluorination 


Reaction 

KJ/gm  of  fuel 

KJ/gm  of  fuel  -t- 
oxidizer 

B(s)  +  3/402(g)  — >  1/28203(1) 

5X 

18 

B(s)  -1-  3/2F2(g)  — >  BF3(g) 

105 

17 

B(s)  -f  l/202(g)  +  l/2F2(g)  — >  OBF(g) 

56 

13 

The  extent  to  which  the  thermodynamic  potential  of  boron/fiuorine  s>sLems  is  realized 
ultimately  depends  on  the  combustion  kinetics.  The  goal  of  the  present  work  is  to  develop  a  kinetic 
model  for  boron  combustion  that  can  be  used  in  helping  to  identify  critical  parameters  and,  in 
conjunction  with  experimental  studies,  to  better  understand  combustion  mechanisms.  Since 
particulate  boron  combustion  is  a  complex  multisiep  proce,ss  involving  both  heterogeneous  and 
homogeneous  chemical  kinetics,  this  is  an  inherently  muupha.sed  task.  Under  an  earlier  project 
sponsored  by  the  Air  Force  Office  of  Scientific  Re.search  models  were  developed  for  the 
homogeneous and  heterogeneous-^'-'’  kinetics  a.ssociated  with  particulate  boron  in  post- 
hydrocarbon  combustion  gases.  The  heterogeneous  processes  treated  included  both  the  high 


icmperaiure  surtacc  buminu  ol  a  relauveiy  clean'  boron  panicle"^  and  chemically  taciliiateu 
gasiiicaiion  ot  ihe  boron  oxide  coaung-*'.  The  overall  goal  of  the  present  Oltice  of  Naval  Research 
effort  is  to  extend  the  model  lor  B/H/O/C  combustion  systems  to  ip.clude  fluorine  chemistrs  The 
specific  objecuves  are  to  characienzc  the  general  mechanistic  behavior  of  these  systems,  to  identify 
the  most  important  gas-phase  reaction  pathways  and.  hence,  to  determine  the  reaction  rate 
parameters  who.se  experimental  and  theoretical  evaluauon  would  significantly  enhance  current 
predictive  capabilities. 

Initial  work  extended  the  earlier  gas-phase  o.xidauon  model  lor  B/O/H/C  combustion 
.systems'"^  to  include  fluorine  chemistry*^.  It  included  a  destnption  of  a  reaction  mechanism  for 
high  temperature  B/H/O/C/F  combustion  systems,  the  re.sulLs  of  kinetic  calculations  illustrating  the 
mechanistic  behavior  of  these  systems,  and  an  analysis  of  the  dominate  reaction  fluxes  and  key 
reaction  pathways.  Here  we  f  ormulate  models  for  the  heterogeneous  gas-surf  ace  kinetics.  The 
study  specif  ically  treats  the  gasification  of  a  liquid  boron  oxide  droplet.  However,  a  preliminary 
summary  of  the  surface  kinetics  model  for  oxidation  of  a  .solid  boron  particle  is  also  reported  for 
completene.ss.  However,  the  latter  model  has  not  been  extensively  studied  and  only  .ser\es  as  the 
starting  point  for  more  detailed  analysis. 

The  remainder  of  this  report  is  organized  as  follows.  Section  2.0  pre.senLs  an  overview  of 
the  single  panicle  combustion  model  and  summarizes  the  governing  equations  for  gasificauon  of  a 
liquid  boron  oxide  droplet.  Section  3.0  reviews  the  gas  phase  oxidation  model.  Section  4.0 
describes  the  heterogeneous  B2O3  gasification  kinetics  model.  Section  5.0  presents  model  results. 
Section  6.0  summarizes  the  work  to  date  and  describes  further  research  efforts.  All  technical 
references  are  listed  in  Section  7.0.  Appendix  A  contains  a  complete  listing  of  the  gas  pha.se 
oxidation  mechanism.  Appendix  B  lists  the  preliminary  heterogeneous  (gas-surface)  reacuon 
mechanism  for  solid  boron  in  H/O/C/F  combustion  environments. 


2.0  SPHERICAL  PARTICLE  COMBUSTION  MODEL 


2. 1  Model  Overview 

The  sphencal  partici'  eombusiion  model  treats  the  combustion  ol  a  sphencal  ( 1- 
dimensional)  particle  m  ..  not  quiescent  oxidizing  environment.  The  key  processes  are  shown 
schematically  in  Figure  1 . 


Liquid  Boron  Oxide  or 
Solid  Boron  Particle 


Figure  1 .  Schematic  ot  Spherical  Particle  Combustion  in  Reactive  Environments. 


The  present  report  summanzes  a  numeneal  model  and  analysis  lor  the  ehemically  laeilitated 
easitication  ot  a  htjuid  boron  oxide  droplet  m  a  high  temperature  quiescent  environment  consisting 
oi  post  combusuon  gases  containing  lluorine.  The  physical  model  is  that  oi  a  liquid  boron  oxide 
droplet  instantaneously  placed  in  a  quiescent  high  temperature  oxidizing  environment  consisung  ol 
equilibnum  produces  ot  lluonnated  derivauve  ot  HMX  and  TMETN  binder/'  Tlie  eltecLs  ol 
temperature,  mixture  composiuon.  and  droplet  diameter  on  the  gasiiication  rate  are  reported.  The 
numerical  model  and  solution  procedure  have  been  described  in  detail  el.sewhere^  and  are  only 
hnetly  desenbed  here.  The  calculations  were  perlormed  under  quasi-steady  conditions  in  which 
the  droplet  diameter  was  held  fixed  in  ume  and  under  transient  conditions  in  which  the  droplet 
diameter  was  allowed  to  regress. 

2  2  Governing  Equations 

The  low  Mach  number  gas-pha.se  tlow  is  de.scribed  as  an  ideal  gas  with  elementary  detailed 
kinetics  and  multi-component  ma.ss  and  thermal  diflu-sion.  The  assumption  ol  isobanc  How 
.ipplies  and  the  momentum  equation  is  .saustied  trivially  by  the  Stelan  ilow  away  Irom  the  panicle 
surtace.  The  a.ssumptions  employed  here  include:  { 1 )  negligible  viscous  dissipauon.  (2)  no 
buoyancy  eficcLs.  (./)  no  torced  convection  on  the  particle,  and  (4)  no  Dulour  eliects.  The  range 
ol  ambient  pre.ssures  to  be  studied  by  the  model  arc  limited  to  tho.se  lor  which  the  ideal  gas 
equation  ol  state  and  kinetics  arc  valid.  Thus,  the  boundary  .separating  the  conden.sed-phase  from 
the  gas-pha.se  will  be  treated  as  a  discrete  interface. 

The  equations  governing  the  gas-phase  arc: 


c)t 


'-dr 


(r^p^Vf.)  =  0 


(1) 


p  1,’  (r  p —  ()  (A— 1 . (j) 

dr 


(2) 


dT, 


r.^~ 


dr 


dT,,  1  d  2.  VI/. 

Px^px^’x  -37  Px^xk^x'^'-p 


dT, 


dr  2  dr  ^  dr  ,  . 
r  k=  I 


xk' 


dr 


G 


+  ^  kikgkf' xk^’’ Xk  =  P  (^) 


k=\ 


Px  = 


pWg 

RTg 


(4) 


In  these  equations,  r  denotes  the  spatial  radial  coordinate:  r.  the  lime:  pi>.  the  mass  density 
of  the  gas-pha.se:  i'^,.  the  radial  velocity  of  the  ga.s-pha.se;  the  mass  fraction  of  the  kih  gas- 


phase  species;  V(>ii,  ihe  dit'tusion  velociiy  oi  w1e  klh  gas-phase  species:  .  die  molar  producuon 

rale  by  chemical  reaction  ot  the  kth  gas-phase  species  per  unit  volume;  W^k.  the  molecular  weight 
ol  the  klh  gas-phase  species:  W,,  the  mean  molecular  weight  ot  the  gas-phase  mixture:  cpif,  the 
constant  pressure  heat  capacity  ot  the  gas-phase  mixture:  cp^i^.  the  constant  pressure  heat  capacity 
ot  the  klh  gas-phase  species:  T^,  the  gas-phase  temperature:  A,c.  the  thermal  conductivity  ol  the 
gas-phase  mixture;  the  specific  enthalpy  of  the  kth  gas-phase  species:  p.  the  pressure;  R.  the 
universal  gas  constant;  and  G.  the  number  of  gas-phase  species. 


The  imerfacial  conditions  at  the  liquid-phase/gas-phase  boundary  {Ig)  for  t  >  0  are: 


r  = 

(5) 

T,=  T, 

(6) 

^l.^k  ~  (^'=l . ,G) 

(7) 

X  (^lak)  = 

k=l 

(8) 

(9) 

In  these  equations,  the  elements  of  (Ok  include  both  the  processes  of  vaporization  and  condensation 
and  of  solution  and  dissolution.  Hence,  includes  both  the  enthalpies  of  solution  and  the  heats  of 
vaporization. 

Tlie  outer  boundary  condition  for  t  >  0  is: 

/-  — >  oo 

(lO) 

lim  y,k  =  y.k^,  (a=i . g) 

r  — >  oo 

(ID 

lim  T,  = 

r  — >  oo 

(12) 

The  initial  conditions  for  the  gas  and  liquid/gas  interface  are: 

<  r  <  oo.  [^.?](=0  ~  (^'=1 . G)  (13) 

()</-<r/^,  [Ti],=()=Tiq  (1 4) 


6 


r  =  Oy, 


[0,?a]/=o  -  (k~\ 


LG) 


115) 


In  these  ecjuaiions,  Ci^k  denotes  the  number  per  unit  surlace  area  ol  the  kth  surtace  complex  at  the 
liquid/gas  interlace. 
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3.0  GAS  PHASE  KINETICS 


This  section  briefly  reviews  the  model  being  used  to  treat  the  gas  phase  kinetics.  More 
complete  discussions  are  given  in  References for  B/H/O/C  systems  and  in  Reference^  for 
B/H/O/C/F  systems. 

3. 1  Gas  Phase  Reactants 

The  reactants  included  in  the  gas  phase  oxidation  model  are  listed  in  Table  2.  They  include 
the  dominant  speciation  in  post-hydrocarbon  combustion  environments  (H/O/C),  boron  oxides  and 
oxyhydrides  that  are  important  intermediates  and  products  in  boron  assisted  hydrocarbon 
combustion  (B/H/O/C),  and  the  additional  boron  fluorides  and  oxyfluorides  that  arise  when 
fluorine  is  included  in  the  system  (B/H/O/C/F). 


Table  2.  B/H/O/C/F  Gas  Phase  Species  List* 


H/O/C 

B/H/O/C 

B/H/O/C/F 

H 

B 

F 

0 

BO 

HF 

OH 

BO2 

BF 

H2O 

B2O2 

BF2 

H2 

B2O3 

BF3 

02 

HBO 

OBF 

CO 

HOBO 

CO2 

HO2 

H2O2 

HCO 

*  HO2,  H2O2  and  HCO  are  used  only  in  describing  H/O/C  chemistry. 


The  reactant  list  in  Table  2  is  not  intended  to  suggest  that  these  will  be  the  only  species  of 
importance  for  an  arbitrary  propellant  formulation.  Rather,  these  reactants  were  selected  based  on 
equilibrium  analyses  of  specific  propellant  formulations  and  a  given  range  of  combustion 
temperature  and  pressure.  The  boron  containing  species  in  Table  2  were  typically  found  to  be 
dominant  (largest  equilibrium  mole  fraction)  for  the  temperatures  of  interest  and  for  a  range  of 
boron  to  oxidizer  mass  ratios.  Notable  species  which  are  not  included  are  the  oxide  B2O.  the  HBO 
i.somer  HOB  and  boron  hydrides  BH  and  BH2.  Additionally,  the  species  H02,  H2O2  and  HCO 
are  treated  only  within  the  kinetic  mechanism  for  H/O/C  systems.  Reactions  between  these  species 
and  boron  species  have  not  been  included.  An  expanded  gas  phase  model  for  B/H/O/C  systems 
which  includes  HCO  reactions,  as  well  as  HOB  and  boron  hydride  has  been  reported  by 
Pasternak.*^ 
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Heals  ot  tormauon  and  entropies  tor  gas  phase  species  are  listed  in  Table  3.  The  heat  ot 
lormaiion  tor  HBO  is  t'rom  Page^.  Similar  values  have  been  reported  by  other  workers,  The 
heal  of  lormaiion  for  HO2  is  from  Hills  and  Howard  •  * .  Thermochemical  constanus  for  the 
remaining  species  are  from  the  J  ANNAF  Tables*^. 


Table  3.  B/O/H/C/F  Thermochemical  Parameters* 


Species 

S298 

g!f 

vmm 

33DDII 

H 

52.10 

27.39 

4.97 

4.97 

4.97 

4.97 

4.97 

0 

59.56 

38.47 

5.23 

5.08 

5.02 

5.00 

4.98 

OH 

9.32 

43.88 

7.15 

7.07 

7.13 

7.33 

7.87 

H2 

.00 

31.21 

6.90 

7.00 

7.07 

7.21 

7.73 

02 

.00 

49.01 

7.01 

7.44 

8.07 

8.35 

8.72 

H2O 

-57.80 

45.10 

8.00 

8.44 

9.22 

9.87 

11.26 

HO2 

3.50 

54.43 

8.36 

9.48 

10.75 

11.37 

12.34 

H2O2 

-32.53 

55.66 

10.41 

12.34 

14.29 

15.21 

16.85 

CO 

-26.42 

47.21 

6.95 

7.14 

7.61 

7.95 

8.41 

CO2 

-94.06 

51.08 

8.91 

10.65 

12.32 

12.99 

13.93 

HCO 

10.40 

53.66 

8.24 

9.28 

10.74 

11.52 

12.56 

B 

133.80 

36.65 

4.97 

4.97 

4.97 

4.97 

4.97 

BO 

.00 

48.60 

7.00 

7.39 

7.86 

8.11 

8.53 

BO2 

-68.00 

54.90 

2.54 

13.05 

13.64 

13.94 

14.43 

B2O2 

-108.99 

57.96 

15.79 

16.80 

18.00 

18.61 

19.65 

B2O3 

-199.80 

67.80 

19.19 

20.52 

22.09 

22.90 

24.26 

HBO 

-60.00 

48.40 

9.38 

10.42 

1 1.66 

12.31 

13.42 

HBO2 

-133.99 

57.27 

11.95 

13.34 

15.01 

15.89 

17.47 

F 

18.86 

37.92 

5.25 

5.19 

5.12 

5.08 

5.02 

BF 

-27.70 

47.89 

7.65 

7.93 

8.28 

8.46 

8.77 

BF2 

-141.(X) 

59.05 

11.73 

12.21 

12.76 

13.03 

13.48 

BF3 

-271.42 

60.77 

15.65 

16.54 

17.59 

18.12 

18.99 

OBF 

-143.99 

53.70 

11.69 

12.34 

13.12 

13.51 

14.18 

HF 

-65.14 

41.51 

6.36 

6.62 

6.98 

7.20 

7.67 

N2 

.00 

45.77 

6.95 

7.08 

7.50 

7.83 

8.32 

*units  are  Kcal/mol  for  AHf,298  and  cal/mol-K  for  S  and  Cp 
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The  heats  ot  formation  in  Table  3  were  used  to  compute  the  gas  phase  bond  energies  listed  in  Table 
4.  These  are  later  used  in  estimating  adsorption  enthalpies  for  the  gas  phase  reactants  in  Table  2  on 
boron  oxide  surfaces. 


Table  4.  Gas  Phase  Bond  Strengths 


Bond 

Bond  Order 

Bond  Strength 

(ev) 

(Kcal/mol) 

F-H 

1 

5.91 

136.2 

F-B 

1 

7.83 

180.5 

F-BF 

1 

5.59 

132.3 

F-BF2 

1 

6.62 

149.4 

F-BO 

1 

7.06 

163.0 

HB=0 

9.23 

225.4 

B=0 

2 

8.37 

193.4 

0BB=0 

8.29 

193.4 

FB=0 

n 

7.62 

175.8 

H0B=0 

2 

7.51 

173.6 

0B0B=0 

2 

7.40 

171.0 

0B=0 

1.5 

5.51 

127.6 

OB-BO 

1 

4.73 

109.0 

HO-BO 

1 

6.22 

143.3 

OBO-BO 

1 

5.75 

131.8 

BOB 

1 

4.80 

110.8 

H-BO 

1 

4.28 

112.1 

H-BH 

1 

4.77 

109.9 

H-B 

1 

3.48 

80.1 

H-O 

1 

4.38 

102.3 

H-OH 

1 

5.12 

119.2 

h-OBO 

1 

5.00 

118.1 

c=o 

2 

11.11 

256.0 

oc=o 

2 

5.54 

127.6 

3.2  Gas  Phase  Reaction  Mechanism 


A  complete  listing  of  gas  phase  reaction  mechanism  is  given  in  Appendix  A.  The  tabulated 
reaction  enthalpies  were  computed  from  reactant  and  product  heats  of  formation  using  the  data 
listed  in  Table  3.  The  reaction  mechanism  is  limited  to  bimolecular  and  termolecular  (three  body) 
reactions.  Higher  order  reactions  such  asA  +  B+  C  =  D  +  E  were  neglected.  Additionally, 
reaction  channels  resulting  in  products  not  listed  in  Table  2,  e.g.  boron  hydrides,  are  excluded  to 
ensure  consistency. 
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The  portion  ot  reaction  mechanism  describing  the  chemistry  involving  H/O/C  species  has 
been  extensive’;-  used  with  good  resulLs  to  describe  the  kineucs  ot  CO/H2/O2  mixtures.'^  Rate 
parameters  tor  the  reactions  involving  boron  containing  species  are  highly  uncenain.  To  date, 
relevant  laboratory  studies  have  been  limited  to  reactions  between  atomic  B  or  BO  with  O2,  CO2 
and  H20'‘^'20  and  to  reactions  between  BF  and  he  oxidants  O2, 0  and  N02.-^  -“  Rate  parameters 
tor  the  remaining  reactions  have  been  estimated.  For  many  of  the  boron  tluoride  and  oxyfluoride 
reactions,  initial  estimates  were  obtained  from  a  propellant  afterburning  study-^.  The  remainder 
were  estimated  using  simple  scaling  relationships  based  on  the  ratios  of  translational,  rotational, 
and  vibrational  partition  functions. 

3.3  Mechanistic  Behavior 


An  example  of  the  homogeneous,  constant  pressure,  adiabatic  kinetics  of  a  B/O/H/C/F 
mixture  is  shown  in  Figure  2.  Only  the  boron-containing  species  are  given.  The  results  show 
OBF  to  be  the  dominant  boron  species.  The  mole  fractions  of  HBO2  and  B2O3  were  typically  less 
than  0.02.  Note  that  the  temperature  rise  is  associated  with  OBF  formation  and  that  all  the  other 
boron-containing  species  appear  to  simultaneously  form  OBF.  The  rate  of  OBF  formation  was 
found  to  be  equal  to  or  faster  than  HBO2  and  B2O3  formation  in  mixtures  without  lluorine.*-- 
These  results  are  typical  of  those  found  for  mixtures  with  an  oxygen/fluorine  mole  ratio  near  unity. 
A  more  extensive  analysis  of  model  results  for  the  gas  phase  chemistry  is  given  in  Reference  6. 


0  1— — I . .  '  .III 

-8  -7  -6  -5 


log  (t/s) 

Figure  2.  Boron  species  and  temperature  profiles  for  an  adiabatic,  constant  pressure  (1  atm) 
homogeneous  system  with  an  initial  temnerature  of  2000  K.  The  initial  mixture  consisted  of 
X(BF)=0.16.  X(BF2)=().1  1,  X(BO)=0.04.  a;E^02)=0.09.  X(OBF)=0.03,  X(CO)=().04, 
X(H2)=0.04.  X(02)=X(H20)=X(HF)=X(N2)=0.  13. 


3.4  Critical  Gas  Phase  Reactions 


The  most  sensitive  reactions  are  reported  in  Table  5.  Most  ot  these  reactions  are  assiKiated 
with  oxidation  ol  initial  reactants  to  OBF.  A  lew  are  associated  with  the  relatively  minor  sequence 
ol  BF  conversion  to  BF2  and  BF3,  which  was  found  to  be  significant  only  lor  O/F  rauos  much 
less  than  one.  The  simultaneous  formation  of  OBF  by  boron  tluorides.  boron  oxides,  and  boron 
oxyhydrides  is  in  conU'ast  to  the  sequential  formation  of  B2O3  found  in  systems  without  tluonne, 
i.e..  BO  — >  BO2  — >  F1B02  — >  B2O3.  It  is  imponant  to  note  that  the  reactions  designated  as 
most  sensitive  in  Table  5  specifically  refers  to  systems  with  both  oxygen  and  lluorine.  When 
tluorine  is  absent  the  critical  reactions  differ  from  those  in  Table  5.  In  this  case  the  earlier  results 
reported  in  References  1-3  are  more  applicable. 


Table  5.  Most  Sensitive  Gas  Phase  Reactions  for  B/H/O/C/F  Systems 


B/H/O/C/F  Reactions 

F  H2  =  HF  -t-  H 

F-I-H20-I- 

OH  -t-  HF 

F  -t-  B2O2 

=  OBF-i-  BO 

F  -H  HBO  = 

OBF  -h  H 

F  -1-  HBO2 

=  OBF  -1-  OH 

BF  0  =  BO  F 

BF  -h  OH  = 

=  OBF  H 

BF  02= 

OBF  +  0 

BF  BO2 

=  OBF  +  BO 

BF2  +  0  = 

OBF  +  F 

BF2  +  H= 

BF-fHF 

BF2  -t-  BO 

=  OBF  BF 

BF2  -t-  HF 

=  BF3  H 

BF2  -I-  BF2  =  BF3  -1-  BF 

BF3  -1-  BO 

=  Bh2  +  OBF 

B/H/O/C  Reactions 

BO  02  = 

BO2  +  0 

BO  -t-  H2  = 

=  HBO  -K  H 

BO  -1-  HBO2  =  B2O3  -1-  H 

B2O2  +  H 

=  BO  -H  HBO 

B2O2  -t-  0 

=  BO  -h  BO2 

B2O2  -1-  OH  =  BO  +  HBO2 

HBO  -1-  OH  =  BO  -1-  H2O 

HBO  -h  OH  =  HBO2  +  H 

H/O/C  Reactions 

H  -H  O2  =  0  -t-  OH 

0-1-  H2  = 

H  +  OH 

4.0  B203(1)  surface  REACTIONS 


A  detailed  descnpiion  ol  the  model  lor  boron  oxide  gasification  in  B/H/O/C  combustion 
systems  has  been  presented  elsewhere  Here,  that  model  is  extended  to  include  reactants 
containing  fluorine. 

4. 1  Global  Reactions 

Global  surface  reactions  involving  gas  phase  reactants  containing  fluorine  were  selected 
using  the  three  criteria  used  in  earlier  work.-*’  The  selection  criteria  are: 

( 1 )  reactions  first  order  in  gas  phase  reactants 

(2)  reactions  whose  reaction  enthalpy  is  less  than  1 10  Kcal/mol 

(3)  reactions  yielding  gas  phase  products  which  are  consistent  with  the  gas  phase  oxidation  model 

Based  on  these  criteria,  four  additional  reactions  were  identified.  They  are  listed  in  Table  6  along 
with  three  reactions  previously  used  in  modeling  B203(l)  gasification  in  B/H/O/C  combustion 
environments. 


Table  6.  Global  Reactions:  B2Q3(1)  +  Z(g)  =  products 


Reaction 

AHr,298 

(Kcal/mol) 

Rl  B 203(1) 0(g)  =  B02(g)  +  B02(g) 

104.2 

R2  B203(1)  -1-  OH(g)  =  B02(g)  +  HOBO(g) 

88.3 

R3  B203(1)  +  H20(g)  =  HOBO(g) -1- HOBO(g) 

89.4 

R4  B203(1) -t- F(g)  =  B02(g) -1- OBF(g) 

68.6 

R5  B203(1)  -(-  F(g)  — >  BO(g)  +  OBF(g)  +  0(g) 

96.2 

R6  B203(1)  -1-  HF(g)  =  HOBO(g)  -i-  OBF(g) 

86.7 

R6  B203(1) BF(g)  =  B202(g) -h  OBF(g) 

74.3 

Reactions  R4  and  R5  are  the  i.soelectronic  analogs  of  reactions  R2  and  R3.  Isoelectronic  analogs 
of  reaction  R6. 
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B203(1)  +  HBO(g)  =  B202(g)  +  HOBO 


AHr  =  1 16.6 


(16) 


B203(l)  +  C0(g)  =  B202(g)  +  C02(g)  AHr=  123.0  (17) 

are  not  included  because  their  reaction  enthalpies  exceeded  1 10  Kcal/mol.  The  model  treats  each  of 
the  reactions  in  Table  6  as  reversible.  The  discussion  that  follows  descnbes  adsorption  and 
desorption  reaction  steps  for  these  reactions  in  the  forward  direction.  The  backward  reactions  are 
treated  by  taking  each  'elementary'  reaction  step  to  be  reversible.  A  consequence  of  this  is  that  the 
model  will  include  both  first-order  (with  respect  to  surface  complexes)  and  second-order 
desorption  processes. 

4.2  Elementary  Reaction  Mechanism 

4.2. 1  First  Order  Adsorption  and  Desorption  Channels 

The  global  reactions  in  Table  6  have  been  represented  by  simple  adsorption  and  desorption 
reaction  steps  using  the  same  model  previously  used  to  treat  B/H/O/C  systems.  This  consists  of 
following  steps. 

Step(l) 

1  he  B2O3  surface  is  taken  to  consist  of  chains  of  BO3  units  and  reactive  surface  sites  are 
represented  by: 


^>B-0-B<q 


(18) 


Step  (2) 

Adsorption  of  a  gas  phase  species  X-Y  is  taken  to  occur  via  bonding  between  the  more 
electronegative  component  (Y)  to  boron.  Depending  on  the  specific  reactant,  the  more 
electropositive  component  (X)  may  bond  to  oxygen. 

X-Y 
1  I 

‘^>B-0-B<q  -I-  XY  — >  ^>B-0— B<q  (19) 

or 

Y-X 

I 

^>B-0-B<q  -t-  XY  — >  ®>B-0— B<q  (20) 

We  also  allow  the  -0-B<  and  X-Y  bonds  to  be  broken  during  adsorption  resulting  in  two 
surface  species: 
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X— Y 


®>B-0— B<q  . >  ®>B-0-X+  %B-Y  (21) 

One  exception  to  the  above  description  is  when  X  is  boron  (e  g.  BF).  In  this  case,  the 
boron  atom  is  taken  to  be  inserted  in  the  >B-0  bond  to  give  >B-B=0.  This  is  the  least 
endothermic  bonding  configuration  and  is  also  consistent  with  B202(g)  as  a  gas  phase 
product  (Reaction  R6). 

Step  (3) 

Based  on  step  (2),  all  adsorption  products  have  the  form  C)>B-Z.  Desorption  is  then  taken 
to  occur  by  breaking  two  B-0  bonds  and  forming  a  B=0  double  bond.  After  extending  the 
BO3  chain  for  clarity,  bond  cleavage  and  formation  during  desorption  consists  of  the 
following: 

surface  species  bond  cleavage  bond  formation 


^>BO 

>B-Z  — > 

®>BO 


^>B.  +  .0-B<o  +  O-B-Z  — >  ^>B-0-B<q  +  0=BZ 


(22) 


For  simplicity,  the  process  will  be  represented  as  ®>B-Z  — >  0=B-Z.  Based  on  steps  (2) 
and  (3),  adsorption  of  gas  phase  reactant  X-Y  yields  two  surface  species.  The  subsequent 
desorption  of  both  complexes  removes  one  B2O3  molecule  from  the  sun'ace. 


The  above  description  is  not  presented  as  a  postulated  mechanism  for  the  surf  ace  reactions 
in  Table  6.  Rather,  given  our  assumptions  concerning  the  B2O3  surlace  structure,  this  model 
gives  a  simple  description  of  bond  formation  and  cleavage  required  to  yield  gas  phase  produces 
which  are  consistent  with  the  global  reactions  in  Table  6.  To  illustrate  this,  the  complete  reaction 
sequence  for  each  reaction  in  Table  6  is  given  in  Table  7  in  terms  of  the  adsorption  and  desorption 
reaction  steps. 
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Table  7.  Model  Reaction  Sequence  tor  B203(l)  Surlace  Reactions 


ID 

Reactants 

Adsorption 

Desorption 

Products 

Products 

O 


1  0>B-0-B<q  +  O  — >  ^>B-0-B<q  — >  ^>B-0.  +  ^>B-0.  — >  0=B0  +  0=B0 

OH 

I 

2  ®>B-0-B<q  +  OH  -->  ®>B-0-B<q  — >  ^>B-OH  +  ^>B-0.  — >  0=B0H  +  0=B0 

H-OH 
I  I 

®>B-0-B<q  +  H20  — >  ®>B-0-B<q  — >  ^>B-OH  +  ®>B-OH  -->  2  0=B0H 

F 
1 

0>B-0-B<q  +  F  — >  ®>B-0-B<q  — >  ^>B-0.  +  ®>B-F  — >  0=BO  +  0=BF 

H-F 
I  I 

^>B-0-B<o  +  HF  — >  ^>B-0-B<q  — >  ®>B-OH  +  ®>B-F  — >  0=B-0H  +  0=BF 

B— -F 

I  1 

®>B-0-B<q  +  BF  — >  ^>B-0-B<q  — >  ®>B-BO  +  ^>B-F  — >  0=B-B=0  +  0=BF 


To  estimate  adsorption  enthalpies,  bond  energies  for  surface  speciation  have  been  estimated 
using  the  gas  phase  bond  energies  in  Table  4.  These  bond  energies  are  listed  in  Table  8  along  with 
the  gas  phase  bond  energies  for  species  with  similar  configurations.  For  D3  we  have  adopted  the 
gas  phase  0-H  bond  strength  in  HOBO.  For  >BB=0  we  have  used  the  B=0  bond  strength  in 
B2O2.  For  the  remaining  complexes,  the  gas  phase  bond  energies  have  been  adjusted  to  refect 
differences  in  boron's  hybridization.  One  difference  between  the  surface  complexes  and  gas  phase 
species  for  D2,  D4.  D5  an  D7  is  the  boron-  oxygen  double  bond  in  the  gas  phase  species  (i.e. 

0=B-Z)  versus  the  two  boron-oxygen  single  bonds  for  the  complex  (i.e.  ^>B-Z).  Consequently, 
in  these  cases  the  bond  energies  for  surface  complexes  have  been  set  equal  to  the  gas  phase  bond 
energies  minus  10  Kcal/mol. 
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Tables.  B203(l)  Surlacc  Species  Bond  Sirengths 


Surface  Complex 
Bond 

Bond  Strength 
(eV) 

Bond  Strength 
(Kcal/mol) 

Gas  Pha.se 

Bond 

Bond  Strength 
(Kcal/mol) 

Di  O>b0-B< 

5.75 

1 10.0 

D2  0>B-F 

6.95 

153.0 

OB-F 

163.0 

D3  0>BO-H 

5.00 

1 18.1 

OBO-H 

1 18.1 

D4  0>B-0H 

6.22 

133.3 

OB-OH 

143.3 

D5  0>B-BO 

4.73 

100.0 

OB-BO 

109.0 

D6  0>BB=0 

8.29 

191.0 

0BB=0 

193.4 

D7  o>b-o. 

5.51 

1 17.0 

OB-O 

127.0 

Adsorption  enthalpies  and  rate  parameters  tor  adsorption  reactions  are  listed  in  Table  9.  For 
the  adsorption  rates  we  have  assumed  a  rate  constant  given  by 


ka  =  kf^fr SoeVRT 


123) 


where 


In  Equation  (23).  Ea  is  the  adsorption  activation  energy  and  sp  is  the  probability  that  a  molecule 
hitting  the  surface  adsorbs.  Rate  parameters  for  adsorption  of  0(g).  OH(g)  and  H20(g)  were 
estimated  previously  and  are  used  here  without  any  change.  Anderson  reports  that  his  preliminary 
re.sulLs  indicate  the  HF  and  H2O  exhibit  about  the  same  reactivity  with  respect  to  B2O3  clusters. 
Consequently.  H2O  adsorption  rate  parameters  have  been  adopted  to  de.scnbe  HF  adsorpuon. 
Similarly,  we  will  also  assume  that  F(g)  has  the  same  rate  parameters  as  OH(g). 
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Table  9.  Firsl-order  Adsorpiion  Rate  Parameters  lor  8203(1)  Surtace 


Adsorpuon  Reacuon 

B203(1)  +  Z(g)  — >  adsorption  produces 

ka  =  koVTs„cl  rtI 

AHf,298 

(Kcal/mol) 

ki) 

tem/s) 

■’0 

Ea 

(Kcal/mol) 

S-1  >B-0-B<  +  0  =  >B-0.  +  >B-0. 

-17 

910. 

0.1 

0.00 

S-2  >B-0-B<  +  OH  =  >B-0.  +  >B-OH 

-33 

882. 

0.1 

0.00 

S-3  >B-0-B<  +  H2O  =  >B-OH  +  >B-OH 

-32 

857. 

0.05 

10.0 

S-4  >B-0-B<  +  F  =  >B-0.  +  >B-F 

-50 

834. 

0.1 

0.(M) 

S-5  >B-0-B<  +  HF=  >B-OH  +  >B-F 

-32 

813. 

0.05 

lO.O 

S-6  >B-0-B<  +  BF=  >B-BO  +  >B-F 

-59 

666. 

0.01 

30.0 

Desorption  enthalpies  are  estimated  using  the  product  heats  of  formation  Table  3  and  the  adsorption 
enthalpies  from  Table  9.  As  an  example,  consider  reactions  R1-R3  in  Table  2.  i.e. 

^>B-0-B<q  +  O  — >  ^>B-0.  +  0>B-0.  — >  B02  +  BO2  (25) 

^>B-0-B<q  +  OH  -->  ®>B-OH  +  ®>B-0.  — >  BO2  +  HOBO  (26) 

^>B-0-B<o+  H2O  — >  0>B-OH  +  ^>B-OH  — >  HOBO  +  HOBO  (27) 

From  (9).  the  desorption  enthalpy  for  >B-0.  is  1/2(H|-Ha(0)l  where  HI  is  the  reaction  enthalpy 
for  the  global  reaction  R1  and  Ha(0)  is  the  adsorption  reaction  for  0(g).  Similarly,  the  desorption 
enthalpy  for  >B-OH  using  (1 1)  is  l/2[H3-Ha(H20)).  As  a  consistency  check,  the  sum  of  the 
desorption  enthalpies  for  >B-0.  and  >B-OH  should  be  clo.se  to  H2-Ha(0H).  A  similar  procedure 
is  u.sed  to  estimate  the  de.sorption  enthalpies  for  >B-F  and  >B-BO. 

First  order  de.sorption  rates  are  taken  to  have  the  form, 

ki  =  Aj  X  exp(-Ei/kT)  (28) 

The  desorption  coefficients  A 1  are  given  approximately  by  the  transition  state  frequency  factor 
kfiT/h.  As  a  first  approximation,  the  desorption  activation  energies  are  taken  to  equal  the 
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dest)rpiion  enthalpies.  A  summarv’  ot  desorption  enthalpies  and  rate  parameters  are  eiven  in  Table 
10.  The  result  ol  these  estimations  is  that  the  desorpuon  rates  lor  the  complexes  >B-Z  are  largely 
independent  oi  Z  whereas.  Irom  Table  S.  adsorption  rates  are  dependent  on  the  gas  phase  reactant. 


Table  10.  First-order  Desorption  Rate  Parameters  tor  B203(l)  Surface 


Desorption  Reaction 

Forward  Rate  Parameters 
ki  =  Aiexpi-Ei/RT) 

S-7  0>B-0  =  0=B-0 

S-8  0>B-OH  =  0=B-OH 
S-9  0>B-F  =  0=B-F 

S-IO  0>B-BO  =  0=B-B0 

AHi  Ai  El 

(Kcal/mol)  ips'*)  (KcaPmol) 

55.6  0.02  xT  55.6 

55.7  0.02  xT  55.7 

56.0  0.02  xT  56.0 

57.3  0.02  xT  57.3 

4.2.2  Reversibility 

It  is  necessary  to  ensure  that  the  reactions  describe  above  are  reversible.  For  the  global 
reactions  in  Table  2  to  be  reversible,  the  adsorption  reactions  in  Table  7  will  also  be  considered 
reversible.  Thus  the  model  also  includes  the  second-order  desorption  reactions  listed  in  Table  1 1. 
Here  again,  we  have  taken  a  simple  Anhenius  rate. 


Table  1 1.  Second-order  Desorption  Rate  Parameters  lor  6203(1)  Surlace 


No. 

Desorption  Reaction 

Forward  Rate  Parameters 

k2  =  A2exp(-E2/RT) 

A2 

E2 

(cm^/complex-s) 

(Kcal/mol) 

S-11 

<->>B-0  -1-  0>B-0  =  B203(1)  -i-  0(g) 

3.0  X  10-5 

5. 

S-12 

0>B-0H  -t-  0>B-0  =  B203(1)  -t-  OH(g) 

3.0  X  10-5 

5. 

S-13 

0>B-OH  -t-0>B-0H  =  B203(1)  -t-  H20(g) 

3.0  X  10-5 

10. 

S-14 

0>B-F  +  0>B-0H  =  B203(1)  -i-  F(g) 

3.0  X  10-5 

5. 

0>B-F  +  0>B-0H  =  B203(1)  -i-  HF(g) 

3.0  X  10-5 

10. 

|S-16 

0>B-F  -1-  0>B-B0  =  B203(1)  +  BF(g) 

3.0  X  10-5 

10. 

From  Table  2.  the  primary  gas  phase  reactants  are  O.  OH.  H2O.  F.  HF  and  BF.  The 
subsequent  gas  pha,se  products  are  BO2,  HOBO.  OBF  and  B2O2.  Thus,  there  are  no  reactions 
between  the.se  products  and  the  B2O3  surface  which  .satisfy  the  three  .selection  c-iteria.  Therefore, 
to  ensure  reversibility,  we  will  assume  that  ad.sorption  of  the.se  products  to  yield  the  onginal 
surface  complex  is  the  only  viable  channel,  i.e.  HOBO  adsorbs  to  yield  a  >B-OH  complex.  B2O2 
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adsorbs  lo  yield  a  >8-60  complex,  etc.  Rale  parameters  lor  these  additional  adsorption  reactions 
where  expressed  in  terms  ot  the  equilibrium  constant  Kp  and  rate  constants  lor  the  adsorption 
reactions  in  Table  9.  The  explicit  lorms  lor  the.se  rates  are  given  in  Table  1 2. 


Table  12.  First-order  Adsorption  Rale  Parameters  lor  8203(1)  Surlace 


No. 

Reaction 

k3 

S-17 

B02(g)  — >  ^>8-0 

K^{B02)\  ^,-11 

S-18 

H802(g)  — >0>8-OH 

jk,.:,K^(H20)RT 

Kj,{HB02)^  ks-n 

S-19 

08F(g)  — >  C)>8-F 

/  k,.^^RT 

Kp{OBF)\  k,.uK^[H20) 

S-20 

B202(g)  ->  0>B-BO 

Kp(BF)  1  k,.j^Ki,[H,OV>T 

k,.;kk,.(,  K^[HF)Kf,{B202)y  k^.n 

5.0  B203(l)  GASIFICATION  MODEL  RESULTS 

5.1  Prototypical  Example 

The  species  proliles  shown  in  Figure  2  are  representative  ot  the  gas-phase  kinetics 
associated  more  with  particulate  boron  combustion  than  with  boron  oxide  gasilication.  In 
particulate  boron  combustion,  boron  oxides  and  fluorides  are  expected  in  the  reacting  mixture.  In 
the  present  study,  the  gas-phase  kinetics  of  interest  are  those  between  the  products  ot  the  liquid 
oxide  surface  reactions  and  the  species  in  the  surrounding  gas-phase  environment.  In  most 
practical  applications  of  boron  as  a  fuel  or  propellant  additive,  the  initial  composition  of  the  gas- 
phase  environment  with  which  the  oxide  coated  boron  particle  reacts  would  be  void  of  boron 
compounds.  Further,  the  major  surface  reaction  products  are  expected  to  be  B2O3  from 
vaporization.  HBO2  from  the  H2O  and  HF  surface  reactions,  and  OBF  from  the  HF  surface 
reaction. 

An  example  of  the  homogeneous  kinetics  in  the  boundary  layer  surrounding  a  gasifying 
boron  oxide  droplet  is  given  in  Figure  3.  Here,  the  composition  of  the  initial  gas-phase 
equilibrium  mixture  was  derived  from  a  JP- 10/air  mixture  with  an  oxidizer/fuel  ratio  ol  0.5  in 
which  I09r  HF  was  substituted  for  5%  O2  and  59c  N2-  To  this  equilibrium  mixture.  39c  of  the 
remaining  N2  was  replaced  by  1%  of  each  of  the  three  surface  products:  B2O3,  HBO2,  and  OBF. 
At  1800  K,  gas-phase  B2O3  is  converted  almost  entirely  to  OBF.  In  addition,  the  formation  ot 
small  amounts  of  HBO2  and  BF3  are  also  predicted.  A  similar  interconversion  of  boron  oxides 
was  reported  previously  in  mixtures  without  fluorine*-^  where  the  equilibrium  shift  reaction  was 
governed  mainly  by  the  overall  reaction  B2O3  +  H2O  =  2HB02.  Also,  note  that  the  time  tor 
formation  of  OBF  and  HBO2  requires  a  few  ms  while  formation  of  BF3  requires  approximately 
one  second.  Thus,  the  time  to  achieve  equilibrium  in  mixtures  typical  of  those  surrounding  a 
gasifying  droplet  can  be  appreciable. 
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Figure  3.  Species  profiles  for  an  adiabatic,  constant  pressure  ( 1  atm)  homogeneous  system 
with  an  initial  temperature  of  1800  K.  The  initial  mixture  consisted  of  X(H)  =  1.22x10'^, 
X(0)  =  3.84x10--*',  X(OH)  =  4.85x10-4  X(H2)  =  1.16x10-5,  X(02)  =  5.02x10*2,  X(H20)  = 
6.80x10-2.  X(H02)  =  6.40x10-2.  X(H202)  =  1.92x10-^,  X(CO)  =  4.4lx I0--\  X(C02)  = 
6.82x10-2,  X(HF)  =  l.(K)xl0->,  X(N2)  =  6.83xl()-l,  X(B203)  =  X(HB02)  =  X(OBF)'’= 

1.00x10-2. 


Species  profiles  for  the  quasi-steady  gasification  of  a  500  iim  diameter  droplet  of  liquid 
boron  oxide  are  reported  in  Figure  4.  The  ambient  temperature  is  1800  K  and  the  pressure  is  1 
atm.  The  initial  gas-phase  mixture  composition  is  the  same  as  that  reported  in  Figure  3  without  the 
substitution  of  the  3%  N2  for  the  three  surface  reaction  products.  In  this  calculation,  the  particle 
diameter  is  fixed  in  time,  and  an  equivalent  amount  of  mass  is  mathematically  supplied  to  the 
particle  as  it  is  consumed.  The  profiles  show  that  OBF  and  HBO2  are  the  dominant  species 
evolving  from  the  particle  surface.  BO2  and  B2O3  are  more  than  2  orders  of  magnitude  smaller  in 
concentration  than  OBF  and  HBO2. 


Figure  4.  Quasi-steady  species  mole  traction  profiles  as  a  function  of  normalized  position 
from  the  surface.  The  ambient  temperature  is  1800  K  and  the  particle  diameter  is  5fX)  )im. 
The  initial  mixture  consisted  of  X(H)  =  1.22x10'^,  X(0)  =  3.84xl()‘^,  X(OH)  =  4.85x10'“^, 
X(H2)  =  1.16x10-5,  X(02)  =  5.02x10-2,  X(H20)  =  6.80x10-2,  X(H02)  =  6.40x10-2, 
X{H202)  =  1.92x10-8,  X(CO)  =  4,41x10-5.  X(C02)  =  6.82x10-2,  x(HF)  =  l.OOxK)-*, 
X(N2)  =  7.13X10-1. 


In  comparison,  the  steady-state  species  profiles  for  a  500  fim  diameter  liquid  boron  oxide 
droplet  in  the  same  gas-phase  environment,  except  with  the  10%  HF  replaced  by  5%  O2  and  5% 
N2.  is  given  in  Figure  5.  The  results  show  HBO2  to  be  the  dominant  gas-phase  boron-containing 
species. 
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Figure  5.  Quasi-steady  species  mole  fraction  profiles  as  a  function  of  normalized  position 
from  the  surface  for  a  gas-phase  environment  without  fluorine.  The  ambient  temperature  is 
1800  K  and  the  diameter  is  500  pm.  The  initial  mixture  consisted  of  X(H)  =  1.22x10'^,  X(0) 
=  3.84x10--^',  X(OH)  =  4.85xiO-‘*.  X(H2)  =  1.16x10'-^.  X(02)  =  l.(K)xl0->,  X(H20)  = 
f).80xl()-“.  X(H02)  =  6.40x10-7.  X(H202)  =  l.92xl()-^.  X(CO)  =  4.41x10--*',  X(C02)  = 
6.82x10-7,  X(N2)  =  7.63x10-'. 


Normalized  sensitivity  c<  efficients  for  the  gasification  rate  with  respect  to  surface  kinetic 
rate  constants  are  reponed  in  Talkie  :3.  This  table  presents  the  sensitivity  gradients,  ln(K)/ln(aj) 
where  K  is  the  gasification  rate  and  aj  is  a  constant  parameter  of  unity  value  that  is  associated  with 
each  of  the  surface  reaction  rate  constants.  The  gasification  rate  is  most  sensitive  to  the  adsorption 
of  HF  and  desorption  of  0>B-OH.  In  contrast,  the  most  sensitive  reaction  in  the  system  without 
fiuorine  is  adsorption  of  H2O  followed  by  desorption  0>B-OH  (see  Table  14).  However,  as 
shown  below,  the  rate  controlling  steps  of  the  gasification  process  are  extremely  dependent  on  the 
size  of  the  droplet,  the  temperature  of  the  environment,  the  total  amount  of  reactants  in  the  mixture 
and  the  fluorine  to  oxygen  ratio  of  the  mixture. 
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Table  13.  Most  Sen.siiive  Surtace  Reaciions  tor  the  Svsiem  of  Fieure  4 


Reaction 

Sensitivity  Coefficient 

B203(1)  +  HF  — >  <->>B-F  +  ^>B-OH 

0.15 

0>B-OH  — >  HB02 

0.14 

0>b-F  ->  OBF 

0.06 

0>B-F  +  0>B-OH  — >  B203(1)  +  HF 

-0.06 

B203(1)  +  H2O  — >  0>B-OH  +  0>B-0H 

0.05 

0>B-OH  +0>B-OH  — >  B203(1)  +  H2O 

-0.04 

Table  14.  Most  Sensitive  Surtace  Reactions 

or  the  System  of  Figure  5 

Reaction 

Sensitivity  Coefficient 

B203(1)  +  H2O  — >  <->>B-OH  +  <->>B-OH 

0.27 

0>B-OH  — >  HBO 

20.19 

0>B-OH  +0>B-OH  — >  B203(1)  +  H2O 

-0.09 

B203(1)  +  OH  — >  0>B-0  +  0>B-OH 

0.02 

The  time  dependent  rate  of  gasification,  in  units  of  surface  area  per  second,  for  the  system 
in  Figure  4  is  shown  in  Figure  6.  In  this  calculation,  the  particle  radius  was  not  fixed  in  time,  but 
allowed  to  regress  as  the  oxide  was  gasified  from  the  surface.  At  21  sec,  the  diameter  of  the 
droplet  is  approximately  100  fim.  In  comparison,  the  time  dependent  gasification  rate  for  the 
system  without  fluorine  is  also  shown  in  Figure  6.  The  results  show  the  gasification  rate  to 
increase  with  addition  of  HF.  Without  fluorine,  the  diameter  of  the  droplet  equals  100  fim  at 
approximately  33  sec.  Note  that  during  the  first  80*^  change  in  diameter,  the  gasification  rate 
drops  by  30-40‘7f .  before  decreasing  significantly  more  during  the  final  stages  of  gasification. 
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Figure  6.  Gasification  rate  as  a  function  of  lime  for  gas  compositions  with  fluorine  (Figure  3) 
and  without  fluorine  (Figure  4)  with  an  initial  temperature  of  1800  K  and  an  initial  particle 
diameter  of  500  um. 


5.2  El'teci  of  Droplet  Diameter 


The  effect  of  the  droplet  diameter  on  the  gasification  rate  is  given  in  Figure  7,  The  range  of 
diameters  studied  was  2(KX)  pm  to  20  pm.  The  gas-composition,  temperature,  and  pressure  of  the 
surrounding  environment  was  the  same  as  that  of  Figure  4.  Gasification  rates  tor  both  fixed 
diameter  droplets  and  for  three  regressing  diameter  droplets  are  presented.  The  symbols  indicate 
the  gasification  rates  for  each  of  the  fixed  diameter  calculations.  The  solid  lines  show  the  time 
dependent  trajectories  of  droplets  with  initial  diameters  of  2000,  1000,  and  5(X)  pm.  fixed 
diameter  results  show  the  gasification  rate  to  increase  rapidly  for  particle  diameters  up  to 
approximately  500  pm  and  to  be  nearly  independent  of  the  diameter  for  droplets  with  diameters 
greater  than  ItKX)  pm.  Classical  theory  predicts  the  gasification  rate  to  be  independent  of  diameter 
if  the  process  is  diffusion  controlled  and  to  be  linearly  dependent  on  the  diameter  if  the  rate  is 
kinetically  controlled.  The  present  results  predict  both  regimes  at  the  smallest  and  largest  diameters 
studied.  However,  the  transition  from  the  kinetically  controlled  to  diffusion  controlled  rate  occurs 
over  an  extremely  wide  range  of  diameters. 


d//um 

Figure  7.  Gasification  rates  of  liquid  B2O3  droplets  as  a  function  of  diameter.  The  iniual 
surrounding  gas-phase  environment  consisted  of  X(H)  =  1.22x10-^’,  X(0)  =  .^.84x10-5, 
X(OH)  =  4.85x1(H,  X(H2)  =  l.l6xU)--T  X(02)  =  5.02x10-2,  X(H20)  =  6.80x10-2, 
X(H02)  =  6.40x10-7,  X(H202)  =  1.92xl()-^,  X(CO)  =  4.41x10--'',  X(C02)  =  6.82x10-2, 
X(HF)  =  l.tXixK)-*,  and  X(N2)  =  7.13x10'*  with  an  initial  temperature  of  1800  K  and 
pressure  of  one  atmosphere.  The  squares  are  the  steady-state  gasification  rates  obtained  when 
the  diameter  is  not  allowed  to  regress.  The  solid  lines  are  the  transient  gasification  rates  as  a 
function  of  diameter  as  the  droplet  is  consumed.  Data  for  the  transient  calculations  are 
pre.sented  for  initial  diameters  of  20()().  KKK).  5(K).  and  200  pm. 
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The  regressing  surface  results  also  reveal  that  tor  initial  droplet  diameters  less  than  about 
5()()  |im.  the  time  dependent  gasilication  rate  closely  follows  the  fixed  diameter  gasification  rates  as 
the  diameter  is  decreased.  However,  for  droplets  with  initial  diameters  greater  than  50()  pm.  the 
rates  are  predicted  to  deviate  substantially  from  the  fixed  diameter  predictions.  This  deviauon  is 
due  to  the  build-up  of  surface  products  in  a  boundary  layer  surrounding  the  particle,  thus 
significantly  altering  the  reactants  attacking  the  surface. 

Surface  reaction  rates  for  the  fixed  diameter  calculations  as  a  f  unction  of  diameter  are 
presented  in  Figure  8.  Desorption  of  the  0>B-OH  surface  complex  is  the  fastest  step  at  small 
diameters  while  adsorption  of  HF  is  the  fastest  step  for  large  diameter  droplets.  For  small  diameter 
droplets,  the  surface  reactions  are  far  from  equilibrium  whereas  for  large  diameters,  the  surface 
reactions  are  nearly  equilibrated. 


Figure  8.  Surface  reaction  flux  profiles  for  the  steady-stale  results  presented  in  Figure  7. 

The  symbols  denote  the  surface  reactions  as  follows;  S-1,  solid  diamond;  S-2,  solid  circle;  S-.^. 
solid  square;  S-4,  inverted  solid  triangle;  S-5.  solid  U'iangle;  S-7.  four  pointed  star;  S-8  live 
pointed  .star;  S-9,  open  pentagon;  S-1 1,  open  diamond;  S-1 2.  open  circle;  S-I3,  open  square; 
S-14.  inverted  open  triangle;  S-L'i.  open  triangle;  S-17,  four  pointed  cross  "x";  S-18.  five 
pointed  cross;  S-iy.  solid  pentagon.  See  Tables  7-10  for  definitions  of  the  surface  reaction 
numbers. 


The  mosi  sensiuve  surlacc  reaciions,  as  a  lunciion  oi  droplet  diameter,  are  reponed  m 
Figure  9.  Ad.sorption  ol  HF  is  the  most  sensitive  reaction  lor  all  diameters  except  20  pm.  At  20 
pm.  ad.sorption  ol  H2O  is  most  .sensitive,  but  decrca.se.s  in  importance  as  the  diameter  is  increased. 
Interestingly,  the  ellect  that  H2O  ad.sorption  has  on  the  gasification  rate  changes  at  a  diameter  ot 
about  150  pm  (as  indicated  by  a  change  in  sign  ol  the  gradient).  For  diameters  le.ss  than  150  pm. 
this  reaction  tends  to  increase  the  gasilication  rale,  whereas  lor  diameters  greater  than  150  pm.  it 
lends  to  decrea.se  the  gasilication  rale.  Desorption  ol  C)>B-OH  follows  a  similar  trend,  enhancing 
gasification  for  diameters  .smaller  than  4(M)  pm  and  inhibiting  gasification  for  diameters  greater  than 
4(K)  pm.  In  contrast,  the  gasification  rate  becomes  increasingly  more  .sensitive  to  the  desorption  of 
the  0>B-F  complex  as  the  diameter  is  increased. 
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Figure  9.  Sensitivity  gradient  profiles,  lnK/ln(aj),  of  the  response  of  the  gasification  rate  to 
variations  in  the  rate  parameters  of  the  surface  reactions.  K  is  the  gasification  rate  and  aj  is  a 
constant  parameter  of  unity  value  that  is  associated  with  each  of  the  surface  reaction  rate 
constants.  The  calculations  are  for  the  steady-state  results  presented  in  Figure  7.  The  symbols 
denote  the  surface  reactions  as  follows:  S-l.  solid  diamond;  S-2.  solid  circle;  S-.l.  .solid  square; 
S-4.  inverted  solid  triangle;  S-5.  solid  triangle;  S-7,  four  pointed  star;  S-8  five  pointed  star;  S- 
9.  open  pentagon;  S-l  1.  open  diamond;  S-l 2.  open  circle:  S-l.^.  open  square;  S-14.  inverted 
open  triangle;  S-15.  open  triangle;  S-17,  four  pointed  cross  "x";  S-18.  live  pointed  cross;  S- 
19.  solid  pentagon.  See  Tables  7-10  for  definitions  of  the  surface  reaction  numbers. 


5.3  Efteci  ot  Gas  Phase  Temperature 


An  increase  in  the  environmental  temperature  resulLs  in  a  change  in  the  composiuon  of  the 
gas  surrounding  the  droplet.  Given  the  same  elemental  composition  as  studied  in  the  previous 
examples,  the  change  in  equilibrium  mixture  composition  with  temperature  is  shown  in  Figure  10. 
The  major  stable  species  HF,  H2O.  CO2,  and  O2  remain  nearly  constant  while  the  troiisient 
radicals.  OH.  O.  and  H,  and  the  dissociation  produces  CO  and  H2  increase  by  1-2  orders  ol 
magnitude  in  concentration  over  the  temperature  range  of  1600  to  2050  K.  Thus,  a  change  in 
temperature  can  alter  the  gasification  rale  by  changing  the  rate  constants  and  by  changing  the 
speciaiion  of  reactants. 


O'  '  ■  '  (p  ‘  'O'  ■  (b-  'Q‘  ■  ^  'O'  ■  -a)  '  ^ 


Figure  10.  Equilibrium  mixture  composition  as  a  function  of  temperature.  The  mixture  is 
based  on  the  environmental  composition  of  Figure  6. 
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The  etleci  of  lemperaiure  on  ihe  sieady-siaie  gasilicaiion  rale  is  shown  in  Fieure  1 1  lor  a  50 
pm  diaiiieier  droplet.  The  caleulalions  were  perlortned  with  a  lixed  diameter  droplet.  The 
easilicaiion  rale  is  predicted  to  increa.se  with  temperature.  The  increase  in  rate  resulLs  Irom  the 
increa.sed  rale  ol  surface  reactions,  the  increa.sed  rate  ol  vapon/.aiion.  and  the  increa.se  in  radical 
concentrations  in  the  surroundinjj  ^tas.  Below  1950  K.  the  jjasilicauon  proce.ss  is  dominated  by 
surtace  reactions.  Between  1950  and  20.50  K.  vaporization  plays  a  equally  important  role  in 
ea.Nilying  the  liquid  boron  oxide  droplet. 


Figure  1 1.  Gasillcaiion  rates  ol  liquid  B2O3  as  a  lunction  of  temperature.  Tlic  initial  mixture 
compositions  are  presented  in  Figure  lo  All  Ihe  calculations  arc  lor  a  lixed  particle  diameter 
ol  50  pm. 


The  rales  ot  ihe  individual  suriace  reactions  as  a  luncuon  ol  temperature  are  presented  in 
Figure  12.  The  desorption  ot  the  0>B-OH  complex  to  torm  HB02(gt  increa.ses  the  most  rapid 
with  temperature  followed  by  the  desorption  of  the  C)>B-F  complex  to  form  OBF(g).  Reaction  ol 
the  oxide  layer  with  HF  and  H2O  are  the  fastest  adsorption  reactions.  How-ever.  surface  reaction 
by  OH  becomes  appreciable  above  1X(X)  K.  The  rales  of  the  second  order  desorption  reacuons  are 
predicted  to  decrease  with  temperature. 
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Figure  12.  Surface  reaction  flux  profiles  for  the  steady-state  results  presented  in  Figure  1 1. 
The  symbols  denote  the  surface  reacuons  as  follows:  S-1.  solid  diamond:  S-2,  solid  circle;  S-3, 
solid  square;  S-4.  inverted  solid  triangle:  S-5,  solid  uriangle;  S-7.  four  pointed  star;  S-8  five 
pointed  star;  S-9.  open  pentagon;  S-1 1.  open  diamond:  S-1 2.  open  circle;  S-1 3,  open  square; 
S-14.  inverted  open  triangle:  S-15.  open  triangle;  S-17,  four  pointed  cross  "x";  S-18.  five 
pointed  cross:  S-iy.  solid  pentagon.  See  Tables  7-11)  for  definitions  of  the  surface  reaction 
numbers 


The  mosl  scnsiuve  surtaee  reaeiions  on  ihe  easilicaiion  process  are  given  in  Figure  13 
The  rapid  increase  in  sensiiiviiy  ol  dcsorpuon  reactions  with  decreasing  lemperaiure  is  indicative  ol 
iheir  rate  limiung  role  below  16(H)  K  tor  the  droplet  and  mixture  composition  under  consideration 
here.  Above  1700  K.  the  rate  limiting  processes  are  shown  to  be  ad.sorpiion  reactions.  However, 
note  that  between  1950  and  2(K)0  K,  an  increase  in  the  rate  ol  the  three  adsorption  prixe-s.ses  ot 
HF.  H2O.  and  OH  changes  trom  one  ol  acceleration  to  one  ot  inhibition.  Tliis  change  in  behavior 
is  accompanied  by  a  .simultaneous  increase  in  the  vaporization  proce.ss. 
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Figure  13.  Sensitivity  gradient  profiles  of  the  response  of  the  gasification  rate  to  variations  in 
the  rate  parameters  of  the  surface  reactions.  The  calculations  arc  for  the  steady-state  results 
presented  in  F'gurc  1 1.  The  symbols  denote  the  surface  reactions  as  follows;  S-1.  solid 
diamond;  S-2.  solid  circle;  S-3.  solid  square:  S-4,  inverted  solid  triangle;  S-5.  solid  triangle;  S- 
7,  four  pointed  star;  S-8  five  pointed  star;  S-9,  open  pentagon;  S-1 1,  open  diamond;  S-12. 
open  circle;  S-1 3,  open  square;  S-14,  inverted  open  triangle;  S-1 5.  open  triangle;  S-1 7.  four 
pointed  cross  "x  ”;  S-18.  five  pointed  cross;  S-19.  solid  pentagon.  Sec  Tables  7-10  for 
definitions  of  the  surface  reaction  numbers 


5.4  ElTecl  ol  Diluuon 


W..  N,  as  .e  on.  n.o.en_n,  " 

varying  ihe  N2  conieni  in  ihe  surrouiKling  m  the  parual 

Figure  14.  A  decrease  in  ihe  N2  conienl  ot  i  e  mi  casUicaiion  rate  as  shown  in 

pressures  ol  each  of  the  reacianis,  which  ^  eiisificaiion  rale  is  increased  by 

Lure  15.  Compared  10  the  pure  vaponzauonra^  O  K^^^^  As  the  nitrogen 

ab;,m  u  .actor  of  5  when  the  No  ,s  toutUy  -P™  ^  ™  ,,er  .surface  reacuons  are 

content  of  the  mixture  is  decreased  trom  10(1^.  .  reduction  in  the  N2  content, 

predicted  to  increase  the  fastest  (see  figure  1 6).  H"*';  J  ^  ,-,,,,.„per  surtace 

L  rates  of  increase  of  fIuT m  s^-  of  HP  to  he  the  nrost 


Fiuure  14  Equilibrium  mixture  composition  as  a 
based  on  the  environmental  composition  ot  Figure  b. 
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Figure  16.  Surface  reaction  tlux  profiles  for  the  steady-state  results  presented  in  Figure  15. 
The  symbols  denote  the  surface  reactions  as  follows:  S-1,  solid  diamond;  S-2.  solid  circle;  S-3, 
solid  square:  S-4,  inverted  solid  uiangle;  S-5,  solid  triangle;  S-7,  four  pointed  star;  S-8  five 
pointed  star;  S-9,  open  pentagon;  S-1 1,  open  diamond;  S-12,  open  circle;  S-1 3,  open  square; 
S-14,  inverted  open  triangle;  S-15,  open  triangle;  S-17,  four  pointed  cross  "x";  S-18.  five 
pointed  cross;  S-19,  solid  pentagon.  See  Tables  7-10  for  definitions  of  the  surface  reaction 
numbers 
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Figure  17.  Sensitivity  gradient  profiles  of  the  response  of  the  gasification  rate  to  variations  in 
the  rate  parameters  of  the  surface  reactions.  The  calculations  are  for  the  steady-state  results 
presented  in  Figure  15.  The  symbols  denote  the  surface  reactions  as  follows;  S-l,  solid 
diamond;  S-2,  solid  circle;  S-3,  solid  square;  S-4,  inverted  solid  triangle;  S-5,  solid  triangle;  S- 
7,  four  pointed  star;  S-8  five  pointed  star;  S-9,  open  pentagon;  S-l  1,  open  diamond;  S-12, 
open  circle;  S-13,  open  square;  S-14,  inverted  open  triangle;  S-15.  open  triangle;  S-17.  four 
pointed  cross  "x";  S-18,  five  pointed  cross;  S-19,  solid  pentagon.  See  Tables  7-10  for 
definitions  of  the  surface  reaction  numbers. 
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5.5  Eft'ect  ot  Fluonne/Oxygen  Mole  Ratio 


The  lluorine/oxN  gen  (F/O)  ratio  in  the  gas-phase  environment  can  be  vaned  in  several 
manners.  In  the  first  two  examples  presented  here,  the  F/O  ratio  was  varied  without  holding  the 
hydrogen  content  fixed.  In  Figure  18.  the  F/O  ratio  of  the  gas-phase  equilibrium  mixture  was 
varied  by  fixing  separately  the  nitrogen,  carbon,  and  hydrogen  atom  balances  and  the  sum  of  the 
oxygen  and  tluorine  atom  balance.  Tlie  mixture  composition  shows  the  CO2  mole  fraction  to 
remain  nearly  the  same  as  the  F/O  ratio  is  increased  until  a  ratio  of  about  1.6.  at  which  point  the  CO 
mole  fraction  increases  rapidly.  The  H2O  and  OH  mole  fractions  decrea.se  with  increasing  F/O  and 
are  essentially  eliminated  from  the  mixture  at  an  F/O  ratio  near  1 .4.  At  the  same  F/O  ratio,  a 
significant  increase  in  the  F-atom  mole  fraction  is  observed.  The  HF  mole  fraction  increases  and 
the  O2  mole  fraction  decrea.ses  with  increasing  F/O  ratio. 


Figure  18.  Equilibrium  mixture  composition  as  a  function  of  iluorine/oxygen  mole  ratio. 
The  mixture  is  based  on  the  environmental  composition  of  Figure  6. 
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The  gusilicaiion  rates  ol  a  50  |im  diameter  droplet  with  the  mixture  compositions  shown  in 
Figure  17  and  an  ambient  temperature  of  18(K)  K  and  pressure  of  1  atm  are  shown  in  Figure  19. 
The  results  show  that  addition  of  a  small  amount  of  fluonne  to  the  mixture  (F/0  =  0.1)  increases 
the  gasification  rate  by  approximately  16%  .  Further  increasing  the  F/0  ratio  to  1.3  increases  the 
gasification  rate  by  23%.  Interestingly,  an  increase  in  the  F/O  from  1.3  to  2.3  increases  the 
gasification  rate  by  78^7r.  The  sudden  increase  in  the  gasification  rate  at  a  F/O  ratio  of  1 .3 
corresponds  to  the  mixture  where  the  H2O  and  OH  concentrations  are  eliminated  from  the  mixture 
and  the  F-atom  concentration  is  significantly  increased. 


Figure  19.  Gasification  rates  of  liquid  B2O3  as  a  function  of  the  F/0  ratio  in  the  surrounding 
mixture.  The  initial  mixture  compositions  are  presented  in  Figure  18.  All  the  calculations  are 
for  a  fixed  diameter  droplet  of  50  |J.m. 
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The  tastest  and  ihe  rale  limiting  surlace  reactions  as  a  lunctiun  ol  F/O  ratio  are  given  in 
Figures  20  and  21.  respectively.  The  reaction  llux  results  show  a  significant  increase  in  the  rate  of 
F-aiom  adsorption  on  the  surtace  at  a  F/O  ratio  of  1.3.  The  sensitivity  gradient  results  show  H2O 
adsorption  to  be  most  sensitive  for  F/O  <  0.25.  HF  adsorption  to  be  most  sensiuve  for  0.25  <  F/O 
<  1.8.  and  F-atom  adsorption  to  be  most  sensitive  for  F/O  >  1.8. 
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Figure  20.  Surface  reaction  flux  profiles  for  the  steady-state  results  presented  in  Figure  18. 
The  symbols  denote  the  surface  reactions  as  follows:  S-1,  solid  diamond:  S-2.  solid  circle;  S-3, 
solid  square;  S-4,  inverted  solid  triangle:  S-5.  solid  triangle;  S-7,  four  pointed  star;  S-8  five 
pointed  star;  S-9,  open  pentagon;  S-1 1,  open  diamond;  S-12,  open  circle;  S-1 3.  open  square: 
S-14.  invened  open  u-iangle;  S-15.  open  triangle;  S-17,  four  ptiinted  cross  "x";  S-18,  five 
pointed  cross;  S-19,  solid  pentagon.  See  Tables  7-10  for  definitions  of  the  surface  reaction 
numbers. 
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Figure  21.  Sensitivity  gradient  profiles  of  the  response  of  the  gasification  rate  to  variations  in 
the  rate  parameters  of  the  surface  reactions.  The  calculations  are  for  the  steady-state  results 
presented  in  Figure  IS.The  symbols  denote  the  surface  reactions  as  follows:  S-1,  solid 
diamond;  S-2,  solid  circle:  S-3,  solid  square:  S-4.  inverted  solid  triangle:  S-5,  solid  triangle:  S- 
7,  four  pointed  star;  S-8  five  pointed  star:  S-9,  open  pentagon;  S-1 1.  open  diamond:  S-12. 
open  circle;  S-1 3,  open  square:  S-14,  inverted  open  triangle:  S-15,  open  triangle;  S-1 7.  four 
fHunted  cross  "x";  S-1 8.  five  ptunted  cross:  S-19.  solid  pentagon.  See  Tables  7-10  for 
definitions  of  the  surface  reaction  numbers. 
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6.0  SUMMARY  AND  CONCLUSIONS 


In  summary,  a  model  has  been  developed  to  describe  the  gasification  ot  a  liquid  boron 
oxide  droplet  in  combustion  environments  having  a  gas  phase  spcciauon  similar  to  that  expected  in 
the  post-combustion  environment  of  advanced  fluoroamino/nuroamino  based  oxidizers.  The 
model  includes  the  gas  phase  B/H/O/C/F  kinetics  formulated  in  previous  work^  and  the  surtace 
kinetics  used  to  treat  oxide  gasification  in  B/H/O/C  environments.^  The  latter  model  has  been 
extended  to  include  heterogeneous  reactions  involving  tluonne  containing  species.  A  senes  of 
calculations  have  been  performed  to  characterize  the  droplets  burning  rate  on  its  size  and  the 
temperature  and  composition  of  the  surround  gas.  The  model  results  suggest  that  tluonne  can  have 
an  important  impact  promoting  oxide  gasification. 

The  addition  ot  fluorine  to  the  gas  phase  mixture  surrounding  a  gasifying  liquid  B2O3 
droplet  at  high  temperatures  has  been  predicted  to  increase  the  rate  of  gasification.  For  the 
conditions  ol  the  present  study,  the  increase  in  the  gasification  rate  ranged  from  40%  to  a  factor  of 
3. 


For  the  conditions  important  to  boron  particle  combustion,  gasification  of  liquid  B2O3  by 
surtace  reactions  was  predicted  to  always  play  an  important  role.  The  most  sensitive  surface 
reaction  was  found  to  be  adsorption  of  HF  except  for  highly  fluorinated  mixtures  (F/0  >  2)  where 
F  atom  adsorption  was  the  most  sensitive  reaction  and  for  mixtures  with  small  quantities  of 
fluorine  (F/0  <  0.2)  where  H2O  adsorption  became  important.  Adsorption  steps  were  found  to  be 
rate  conu'olling  tor  all  conditions  studied  except  at  the  lowest  temperature  ( 1 6(X)  K)  where 
desorption  became  the  rate  limiting  process.  Thus,  the  assumptions  of  considering  only  reactions 
first  order  in  gas  pasha  reactants  was  found  to  be  appropriate  for  this  study. 


The  results  also  indicate  that  gasification  will  be  kinetically  controlled  for  particle  diameters 
relevant  to  boron  combustion.  Transition  of  the  gasification  rate  from  a  purely  kinetically 
controlled  process  to  one  controlled  totally  by  diffusion  was  found  to  occur  over  a  wide  range  of 
diameters.  Also,  boundary  layer  effects  were  found  to  be  significant  for  an  initial  panicle  diameter 
greater  than  5(K)nm.  The  implications  of  these  results,  particularly  with  regard  to  their  transient 
behavior,  are  important  to  the  interpretation  of  particle  experimental  data. 

Recently,  experimental  studies  of  the  interaction  between  HF  and  small  boron  oxide  cluster 
ions  have  become  available. One  interesting  result  of  this  data  is  the  suggestion  that  formation  of 
an  BHOF  surface  complex  may  be  an  important  adsorption  channel.  The  significance  on  this  data 
is  currently  being  investigated  and  will  be  subsequently  be  incorporated  into  the  model. 
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APPENDIX  A.  GAS  PHASE  REACTION  MECHANISM 

Table  15.  B/O/H/C/F  Reacuen  Mechanism 
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Table  16.  B/O/H/C  Reacuon  Mechanism 


No. 

Reacuon 

A 

~  1 

G-3y 

B+02=B0+0 

0.717  X  10'-^ 

0.0 

G-40 

B+0+M=B0+M 

0.1 10  X  10*^ 

0.0 

HrH 

G-41 

BO+02=B02+0 

0.421  X  10>3 

0.0 

G-42 

B0+0+M=B02+M 

0.109  X  10>^ 

0.0 

-1987.0 

G-43 

B+B02=2B0 

0.361  X  K)!-* 

0.0 

0.0 

G-44 

B0+B02+M=B203+M 

0.181  X  10>4 

0.0 

-1987.0 

G-45 

2B02=B203+0 

0.602  X  IO*J 

0.0 

9935.0 

G-46 

B+OH=BO+H 

0.602  X  10*4 

0.0 

0.0 

G-47 

B0+0H=B02+H 

0.241  X  10*^ 

0.0 

0.0 

G-48 

B0+0H+M=HB02+M 

0.363  X  10*5 

0.0 

-1987.0 

G-49 

B0+H20=HB02+H 

0.602  X  10'* 

0.0 

9935.0 

G-5() 

B02+0H=HB02+0 

0.181  X  10*3 

0.0 

993.5 

G-51 

B02+H+M=HB02+M 

0.181  X  lO'^* 

0.0 

-1987.0 

G-52 

B02+H2=HB02+H 

0.181  X  10*3 

0.0 

2981.0 

G-53 

HB02+0H=B02+H20 

0.120  X  10*3 

0.0 

1987.0 

G-54 

B203+H20=2HB02 

0.602  X  10*^ 

0.0 

1 1922.0 

G-55 

B0+H+M=HB0+M 

0.109  X  10*6 

0.0 

-1987.0 

G-56 

B0+H2=HB0+H 

0.452  X  102 

3.53 

3160.0 

G-57 

B0+0H=HB0+0 

0.164  X  104 

2.76 

5015.0 

G-58 

HB0+0=B02+H 

0.482  X  10*4 

0.0 

0.0 

G-59 

HB0+0H=B0+H20 

0.482  X  10*4 

0.0 

0.0 

G-60 

HB0+0H=HB02+H 

0.482  X  10*4 

0.0 

0.0 

G-61 

HB0+0H=B02+H2 

0.602  X  104 

0.0 

69942.0 

G-62 

HB0+02=B02+0H 

0.602  X  104 

0.0 

69942.0 

G-63 

HB0+0+M=HB02+M 

0.363  X  1021 

-0.5 

50072.0 

G-64 

2B0+M=B202+M 

0.363  X  10*4 

0.0 

-1987.0 

G-65 

B202+H=B0+HB0 

0.602  X  10*3 

0.0 

0.0 

G-66 

B202+0=B0+B02 

0.361  X  10*4 

0.0 

0.0 

G-67 

B202+0H=B0+HB02 

0.361  X  10*4 

0.0 

0.0 

G-68 

B202+0H=B02+HB0 

0.602  X  104 

0.0 

69942.0 

G-69 

B202+02=2B02 

0.602  X  104 

0.0 

80076.0 

G-70 

B02+C0=B0+C02 

0.301  X  10*4 

0.0 

1987.0 

G-71 

B+H20=HB0+H 

0.241  X  10*5 

0.0 

2682.0 

G-72 

B+C02=B0+C0 

0.422  X  10** 

0.0 

0.0 

G-73 

HB0+B02=HB02+B0 

0.182  X  !0*3 

0.0 

990.0 

G-74 

B02+HB0=B203+H 

0.182  X  10*3 

0.0 

990.0 

G-75 

B0+HB02=B203+H 

0.482  X  10*3 

0.0 

0.0 
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Table  17,  O/H/C  Reacuon  Mechanism 


No. 

Reacuon 

A 

r>f 

Ea.t 

G-76 

H+02=0+0H 

0.192  X  1()>5 

0.0 

16440.0 

G-77 

0+H2=H+0H 

0.508  X  10-^ 

2,67 

6292.0 

G-78 

0H+H2=H+H20 

0.216  X  10^* 

1.51 

3430.0 

G-79 

20H=0+H20 

0.123  X  lO-'i 

2.62 

-1878.0 

G-8() 

H2+M=2H+M 

0.457  X  1020 

-1.4 

104400.0 

G-81 

20+M=02+M 

0.617  X  10>^ 

-0.5 

0.0 

G-82 

0+H+M=0H+M 

0.472  X  10*^ 

-1.0 

0.0 

G-83 

H+0H+M=H20+M 

0.225  X  102-'^ 

-2.0 

0.0 

G-84 

H+02+M=H02+M 

0.617  X  1020 

-1.42 

0.0 

G-85 

H02+H=H2+02 

0.663  X  10>*» 

0.0 

2126.0 

G-86 

H02+H=20H 

0.169  X  1015 

0.0 

874.0 

G-87 

H02+0=0H+02 

0.181  X  lO*-^ 

0.0 

-397.0 

G-88 

H02+0H=H20+02 

0.145  X  lOG 

-1.0 

0.0 

G-89 

2H02=H202+02 

0.302  X  10>^ 

0.0 

1390.0 

G-9() 

H202+M=20H+M 

0.120  X  10>8 

0.0 

45500.0 

G-91 

H202+H=H20+0H 

O.KK)x  lO'-* 

0.0 

3590.0 

G-92 

H202+H=H2+H02 

0.482  X  10>4 

0.0 

7948.0 

G-93 

H202+0=0H+H02 

0.955  X  107 

2.0 

3970.0 

G-94 

H202+0H=H20+H02 

0.700  X  10‘^ 

0.0 

1430.0 

G-95 

C0+0+M=C02+M 

0.251  X  1014 

0.0 

-4541.0 

G-96 

C0+02=C02+0 

0.253  X  U)!-^ 

0.0 

47690.0 

G-97 

C0+0H=C02+H 

0.150X  108 

1.3 

-765.0 

G-98 

C0+H02=C02+0H 

0.602  X  10>4 

0.0 

22950.0 

G-99 

HCO+M=H+CO+M 

0.186  X  10'» 

-1.0 

17000.0 

G-KM) 

HC0+02=C0+H02 

0.758  X  1013 

0.0 

410.0 

G-lOl 

HC0+H=C0+H2 

0.723  X  10'4 

0.0 

0.0 

G-l()2 

HC0+0=C0+0H 

0.302  X  10*4 

0.0 

0.0 

G-l()3 

HC0+0H=C0+H20 

0.302  X  1014 

0.0 

0.0 
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APPENDIX  B.  B(s)  SURFACE  OXIDATION  MECHANISM 


B-1.  Global  Reactions 

Global  reactions  tor  paruculate  boron  surtace  oxidation  in  B/H/O/C/F  combustion  systems 
are  listed  in  Table  1 8.  These  reactions  were  selected  using  three  criteria:  ( 1 )  reactions  tirst  order  in 
gas  phase  reactants  (2)  reactions  whose  reaction  enthalpy  is  less  than  150  Kcal/mol  (3)  reactions 
yielding  gas  phase  products  which  are  consistent  with  the  B/H/O/C/F  gas  phase  oxidation  model. 
Reactions  R  i  -R 1 5  are  the  surtace  reactions  previously  used  to  u-eat  B/H/O/C  systems.  Reactions 
RI6-R2O  are  the  new  surtace  reactions  which  arise  when  fluorine  is  added  to  the  system. 


Table  18.  Global  Reactions  B(s) -f  Z(g)  =  Products 


Reaction 

AHr(Kcal/mol) 

R]  B(s)-t-0(g)  =  BO(g) 

-59.6 

R2  B(s)  +  02(g)  =  B02(g) 

-68.0 

R3  B(s)  +  02(g)  =  BO(g)  +  0(g) 

59.6 

R4  2B(s) -i- 02(g)  =  B202(g) 

-109.0 

R5  2B(s) -1- 02(g)  =  2B0(g) 

0.0 

R6  B(s)  -h  OH(g)  =  HBO(g) 

-56.7 

R7  B(s)  +  OH(g)  =  BO(g)  -H  H(g) 

42.8 

R8  B(s)  -I-  H20(g)  =  BO(g)  +  H2(g) 

57.8 

R9  B(s)  -1-  H20(g)  =  HBO(g)  +  H(g) 

.5 

RlO  B(s)  -1-  B02(g)  =  B202(g) 

-41.0 

Rll  B(s)  +  B02(g)  =  2B0(g) 

68.0 

R12  B(s) -t- B203(g)  =  B202(g)  +  BO(g) 

90.8 

R13  B(s) -t- HOBO(g)  =  BO(g) -1- HBO(g) 

86.6 

R14  B(s)  +  HOBO(g)  =  B202(g)  +  H(g) 

77.1 

R15  B(s) -I- C02(g)  =  BO(g)  +  CO(g) 

67.7 

R16  B(s)-i-F(g)  =  BF(g) 

-46.7 

R17  B(s)-i-HF(g)  =  BF(g)-HH(g) 

89.5 

R18  B(.s) -1- BF2(g)  =  BF(g)  +  BF(g) 

85.6 

R19  B(si -t- BF3(g)  =  BF(g) -1- BF2(g) 

102.7 

R20  B(«)-^OBF(g)  =  BO(g) -1- BF(g) 

116.3 

B-2.  Adsorption  and  Desorption  Channels 

Surtace  speciation  is  listed  in  Table  20  and  first-order  adsorption  rate  paramters  are  given  in 
Table  21 .  The  adsorption  enthalpies  are  expresses  in  terms  of  the  heats  of  formation  given  in  Table 
20  for  the  surface  complexes.  First-order  desorption  channels  and  rate  parameters  are  listed  in 
Table  22. 
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Table  21.  First-order  Adsorption  Rate  Parameters  for  Boron  Surface* 


Adsoprtion  Reaction 

AHa 

Ka 

^>0 

Ea 

(Kcal/mol) 

(cm/s) 

(cal/mol) 

H  = 
O  = 
^02(g 
^02(g 
02(g) 
OH(g 
OH(g 
-H2O1 
B02( 

-I-  BO2 

+  B2C 
+  HOI 
-t-HOI 
+  CO2 
-hCO( 
-hBO( 
+  B2C 
+  HB( 
+  F(g) 
+  HFf 


ST 

m 

ST 


c) 

c)  4-  HBO(g) 
2(c) 

(0  +  BO(g) 
02(g)  +  BO(g 
BO(c)  +  BO(g 


2- 52.1 
1-59.6 
09.0 

1 

l-t-59.6 

1- 42.8 

3- 0.3 

2- 2.2 
54-68.0 
14-68.0 
18 

24-134.0 

54-186.1 

14-67.7 

44-26.5 

1 

54-109.0 

34-6O.O 

6 

64-52.) 

6-^H2 

6-27.7 

6-141.0 

6 

6 


everse  ol 

everse  of 

everse  of 

182 

157 

157 

'00 

>00 

149 


BO(c)  desorpuon 
B202(c)  desorption 
HBO(c)  desorption 


[).02 

f).0062 

0.0062 

0.004 

0.004 

0.0025 


1000.0 

10,000.0 

10,000.0 

3000.0 

3000.0 

4000.0 


everse  of  BF(c)  desorption 
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*  Hi  represents  the  heal  ol  t'ormaiion  lor  surface  species  Cj 
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Table  22.  B(s)  Surface  -  Desorption  Rale  Parameters* 


Complex 

Desorption  Channels 

AHi 

(Kcal/mol) 

Ai 

(ps-1) 

El 

(Kcal/mol) 

BH(c) 

B(s)  +  H(g) 

52.1-H2 

reverse  S- 1 

BO(c) 

BO(g) 

B(s)  +  0(g) 

-HI 

59.6-Hl 

().{)2xT 

reverse  S-2 

-HI 

B202 

B202(g) 

2B(s)  +  02(g) 

-I()9.0-H5 

-68.O-H5 

().02xT 

reverse  S-3 

-1()9.{)-H5 

HBO(c) 

HBO(g) 

BO(g)  +  H(g) 
B(s)  +  OH(g) 

-6O.O-H3 

52.I-H3 

9.3-H3 

().02xT 
0.02xT 
reverse  S-7 

-6O.O-H3 

52.I-H3 

BCO(c) 

B(s)  +  CO(g) 

-26.4-H4 

reverse  S- 1 5 

BF(c) 

BF(g) 

B(s)  +  F(g) 

-H6 

18.97-H6 

().02xT 

reverse  S-19 

-H6  j 

*  Hi  represents  the  heat  of  formation 

for  surface  species  Ci 
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